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Abstract 

For the generalized master equations derived by Karrlein and Grabert for the 
microscopic model of a damped harmonic oscillator, the conditions for purity of 
states are written, in particular for different initial conditions and different types 
of damping, including Ohmic, Drude and weak coupling cases, Agarwal and 
Weidlich-Haake models. It is shown that the states which remain pure are the 
squeezed states with constant in time variances. For pure states, the generalized 
nonlinear Schrodinger-type equations corresponding to these master equations 
are also obtained. Then the condition for purity of states of a damped harmonic 
oscillator is considered in the framework of Lindblad theory for open quantum 
systems. For a special choice of the environment coefficients, the correlated 
coherent states with constant variances and covariance are shown to be the only 
states which remain pure all the time during the evolution of the considered 
system. In Karrlein-Grabert and Lindblad models, as well as in the considered 
particular models, the expressions of the rate of entropy production is written 
and it is shown that the states which preserve their purity in time are also the 
states which minimize the entropy production and, therefore, they are the most 
stable ones under evolution in the presence of the environment and play an 
important role in the description of decoherence phenomenon. 
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1 Introduction 

In the last two decades, more and more interest has arisen about the search for a 
consistent description of open quantum systems [1-5] (for a recent review see Ref. 
[6]). Dissipation in an open system results from microscopic reversible interactions 



between the observable system and the environment. Because dissipative processes 
imply irreversibility and, therefore, a preferred direction in time, it is generally thought 
that quantum dynamical semigroups are the basic tools to introduce dissipation in 
quantum mechanics. In Markov approximation and for weakly damped systems, the 
most general form of the generators of such semigroups was given by Lindblad [7]. 
This formalism has been studied for the case of damped harmonic oscillators [6, 8-12] 
and applied to various physical phenomena, for instance, to the damping of collective 
modes in deep inelastic collisions in nuclear physics [13]. A phase space representation 
for open quantum systems within the Lindblad theory was given in [14, 15]. Important 
progress beyond the limitations of the weak coupling approach was made in describing 
quantum dissipative systems by using path integral techniques [16, 17]. 

In the present study we are also concerned with the observable system of a har- 
monic oscillator which interacts with the environment. We discuss under what con- 
ditions the open system can be described by a quantum mechanical pure state. In 
Sec. 2 we present the generalized uncertainty relations and the correlated coherent 
states, first introduced in [18], which minimize these relations. In Sec. 3 we consider 
generalized quantum master equations derived by Karrlein and Grabert in Ref. [17] 
for the microscopic model of a harmonic oscillator coupled to a harmonic bath [19] by 
using the path integral and we obtain conditions for the purity of states, in partic- 
ular for different initial conditions and different types of damping, including Ohmic, 
Drude and weak coupling cases, Agarwal and Weidlich-Haake models. We show that 
the states which satisfy the conditions of purity are the pure squeezed states with 
well-determined constant in time variances. For pure states, we also derive the gener- 
alized Schrodinger-type nonlinear equations corresponding to these master equations. 
The Lindblad theory for open quantum systems is considered in Sec. 4. For the one- 
dimensional harmonic oscillator as an open system, we show that for a special choice of 
the diffusion coefficients, the correlated coherent states, taken as initial states, remain 
pure for all time during the evolution. In some other simple models of the damped 
harmonic oscillator in the framework of quantum statistical theory [20, 21], it was 
shown that the Glauber coherent states remain pure during the evolution and in all 
other cases the oscillator immediately evolves into mixtures. In this respect we gen- 
eralize this result and also our previous result from [10] as well as the results of other 
authors [22], obtained by using different methods. In Sec. 5 we introduce the linear 
entropy, we present its role for the description of the decoherence phenomenon and we 
also derive the expressions of the rate of entropy production. We show in Lindblad 
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and Karrlein-Grabert models that the correlated coherent states, respectively the pure 
squeezed states, which fulfill the condition for purity of states, are also the most stable 
states under evolution in the presence of the environment and make the connection 
with the work done in this field by other authors [23-28]. Finally, a summary is given 
in Sec. 6. 



2 Generalized uncertainty relations 

In the following we denote by oaa =< A 2 > — < A > 2 the dispersion of the operator 
A, where < A >= a a = Tr(/L4),Trp = 1 and p is the statistical operator (density 
matrix). By a AB = (1/2) < AB + BA > — < A X B > we denote the correlation 
(covariance) of the operators A and B. Schrodinger [29] and Robertson [30] proved 
that for any Hermitian operators A and B and for pure quantum states the following 
generalized uncertainty relation holds: 

For the particular case of the operators of momentum p and coordinate q the uncer- 
tainty relation (1) becomes 

a = a pp a qq -a 2 pq >—. (2) 

This result was generalized for arbitrary operators (in general non-Hermitian) and for 
the most general case of mixed states in [18]. The inequality (2) can also be represented 
in the form: 

h 2 

OppOqq > 4 ^ _ (3) 

where 

r = apq (4) 

\J G pp a qq 

is the correlation coefficient. The equality in the relation (2) is realized for a special 
class of pure states, called correlated coherent states or squeezed coherent states, which 
are represented by Gaussian wave packets in the coordinate representation. These 
minimizing states, which generalize the Glauber coherent states, are eigenstates of an 
operator of the form: 

1 IT Tj 
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with real parameters r and r], \r\ < l,r) — y/cr qq . Their normalized eigenfunctions, 
denoted as correlated coherent states, have the form: 

with a a complex number. If we set r = and r\ = (h/2mu) 1 ^ 2 , where m and to are the 
mass and, respectively, the frequency of the harmonic oscillator, the states (6) become 
the usual Glauber coherent states. In Wigner representation, the states (6) look: 

1 (q — a ) 2 2rf 2r 

W r , v (p, q) = — exp[- ^ 2(1 _ 9 r2) - -s-(p - * P ) 2 + %{l _ r2)l/2 {q - * q )(P ~ <r P )], (7) 

where a q and a p are the expectation values of coordinate and momentum, respectively. 
This is the classical normal distribution with the dispersion 

2 ^ ^ r 

ar qq = V , °VV = 4^2(]_ _ f .2^ ' °"P9 = 2(1 - r 2)l/2 ^ 

and the correlation coefficient r. The Gaussian distribution (7) is the only positive 
Wigner distribution for a pure state [31]. All other Wigner functions that describe 
pure states necessarily take on negative values for some values of p, q. 

In the case of relation (1) the equality is generally obtained only for pure states. 
For any density matrix in the coordinate representation (normalized to unity) the 
following relation must be fulfilled: 

1 



v 



Trp 2 < 1. (9) 



The quantity v characterizes the degree of purity of the state: for pure states v — 1 
and for mixed states v > 1. In the language of the Wigner function the condition (9) 
has the form: 

i = 2jrh J W 2 {p, q)dpdq < 1. (10) 

Let us consider the most general mixed squeezed states described by the Wigner 
function of the generic Gaussian form with five real parameters: 

W<J> ' ^ = 2^/^ ex v{-^ a pp(Q ~ a i) 2 + a ««(P - a p f ~ 2o pM - a g)(P ~ a p)}}' ( U ) 

where a is the determinant of the dispersion (correlation) matrix ( ° w ° vq j . Here, 

\ Cpg &qq / 

a is also the Wigner function area - a measure of the phase space area in which 
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the Gaussian density matrix is localized. For Gaussian states of the form (11) the 
coefficient of purity v is given by 

v = ~Jd. (12) 

The inequality (2) must be fulfilled actually for any states, not only Gaussian. Any 
Gaussian pure state minimizes the relation (2). For a > % 2 /A the function (11) corre- 
sponds to mixed quantum states, while in the case of the equality a = h 2 /4 it takes 
the form (7) corresponding to pure correlated coherent states. 

The degree of the purity of a state can also be characterized by the quantum 
entropy (we put the Boltzmann's constant ks — 1): 

S= -Tr(plnp) = - < hip > . (13) 

For quantum pure states the entropy is identically equal to zero. It was shown [12, 32] 
that for Gaussian states with the Wigner functions (11) the entropy can be expressed 
through a only: 

_ u + 1 u+1 u-1 v-\ 2 

S = —— In— —In——, u=-y/a. (14) 

2 2 2 2ft v 1 

3 Generalized quantum master equations 

In the framework of the standard microscopic model [19, 32-34] for the damped har- 
monic oscillator, it was shown in Ref. [17] that in general there exists no exact master 
equation for the damped harmonic oscillator 

= Cp{t) (15) 

with a dissipative Liouville operator C describing the dynamics of the oscillator and 
independent of the initial states. For specific initial preparations the time evolution 
is described exactly by a time- dependent Liouville operator and the resulting master 
equation for the damped harmonic oscillator with the Hamiltonian 

1 2 Mu 2 2 
Ho = 2M P + — q ' (16) 

corresponding to this Liouville operator (given by Eq. (40) of Ref. [17]), has the 

following general form ({, } denotes the anti-commutator of two operators): 

% M M 2 

-2p>(f)[?,W(*)}] + tfD q (t)\p, [q,p(t)}\ ~ jfD p {t)[q, [q, p(t)}\. (17) 
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The dependence on uj is included in the coefficients of commutators. 

A) For the so-called thermal initial condition [36] , which can be used to describe 
initial states of the entire system (oscillator and bath) resulting from position mea- 
surements, D q (t) and D p {t) can be written as 

D g (t) = lq {t) <g 2 > D p (t) = l P (t)^^- (18) 

Here j q (t),^ p (t) and the equilibrium variances < q 2 > and < p 2 > are given in terms 
of the equilibrium coordinate autocorrelation function < q(t)q > . 

We now derive the necessary and sufficient condition for p(t) to be a pure state 
for all times. p(t) is a pure state if and only if Trp 2 (t) = 1. This is equivalent with 
(d/dt)Trp 2 (t) = for all times, i. e. with the condition Tr(p(t)£p(t)) = 0. With the 
explicit form of Cp(t) given by the right-hand side of Eq. (17) and using the relations 
p 2 (t) = p(t) and p(t)Ap(t) = Tr (p(t) A) p(t), we obtain the following condition for a 
state to be pure for all times : 

h 2 

M 2 D p {t)a qq {t) - MD q (t)<j pg {t) - — 7p (t) = (19) 
and by inserting the expressions (18): 

7p (t) < p 2 > a qq {t) - [M lq {t) <q 2 > -^±}<j pq (t) - ^ lp (t) = 0. (20) 

B) For factorizing initial conditions, namely if the initial density matrix of the 
entire system factorizes in the density matrix of the oscillator and the canonical density 
matrix of the unperturbed heat bath [34, 36-38], i. e. if the oscillator and bath are 
uncorrelated in the initial state, then the resulting master equation is equivalent to the 
result by Haake and Reibold [37], who derived it directly from microscopic dynamics 
and by Hu, Paz and Zhang [39] from the path integral representation. The condition 
for purity of states for these master equations has also the form (20), where now the 
coefficients are given by Eqs. (87), (89) in Ref. [17]. 

For a pure state p(t) = \ip{t) >< ip{t)\, we can obtain from Eq. (17) the evolution 
equation for the wave function ip(t) as an equation of the Schrodinger-type 

(21) 

Taking into account the condition for purity of states (19), we find the non-Hermitian 
Hamiltonian 

v 2 M 1 

H ' = 2M + Y lq{t)q2 + 2 lp{t){qP + a ' pit)q ~ aq{t)p) 

i M i M 2 

+^D q (t)(p - a p (t))(q - a q (t)) - —D p (t)(q - a q (t)) 2 , (22) 
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which depends on the wave function ip{t) via the expectation values a q and a p , i. e. 
this Schrodinger-type equation is nonlinear. 

The master equations considered up to now in this Section are exact. We now 
consider particular types of damping for which the dynamics can be described in 
terms of approximate Liouville evolution operators, valid for arbitrary initial states 
[17]. Then the evolution operator is time-independent and the master equation for the 
density matrix obeys Eq. (17), where D q (i) = D q and D p (t) = D p read 

D q = %<<! 2 >-^, A> = 7 P ^, (23) 

with time-independent coefficients 7 g and 7 P [17]. 

1) For strictly Ohmic damping, 7 P = 7 is the Laplace transform of the damping 
kernel [17, 19] of the model and 7 ? = Uq. In this case we do not have a well-defined 
Liouville operator since < p 2 > and, therefore, the coefficients D q and D, p given by 
(23) are logarithmically divergent. The condition for purity of states for this strictly 
Ohmic damping is similar to the relation (20), only now the coefficients are constant: 

7 < P 2 > <r„(t) - (Muj 2 <q 2 > -^±)a pq (t) - 1 7 = 0. (24) 

2) A more realistic case is the so-called Drude damping. For a sufficiently large 
Drude parameter ujd and sufficiently high temperature ksT ^> Try, the oscillator dy- 
namics can be described by an approximate Liouville operator with the coefficients 

7 9 = a 2 + rf, 7 P = 2a, (25) 

where a and r\ depend on 7, u and ujd and then the condition for purity of states is 

< v 2 > h 2 

2a<p 2 > a qq (t) - [M(a 2 + rf) < q 2 > j^K q (t) ~ y« = 0. (26) 

For a pure state, the Schrodinger equation (21) corresponding to the master equation 
with Ohmic damping has the Hamiltonian (22), with the coefficients given by (23) and 
with 7 P = 7, 7 g = u>q. A similar result holds for the Drude damping, when we take 
the coefficients (25). 

3) In Ref. [17] it is shown that in the limit of weak damping the general master 
equation has the following form: 

p(t) = + Y^ 2 ° + ^^(O] 

-^[?, W(*)}] " jj^\P, hP(t)}} - |% [q,p(t)}}- (27) 
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This equation is given in terms of four dissipation coefficients: 7 S leads to a frequency 
shift and may be absorbed by renormalizing cu , 7 C is the classical damping coefficient 
and the coefficients K s and K c depend on the temperature. K s can be calculated 
analytically only in certain cases. One of these is the Drude model. The master 
equation (27) is a generalization of the Agarwal equation [40]: 

= -^"[^ + ^q\ P (t)} - -[q, {p,p(t)}] - K ^coth(^)[g, [g, p(t)]], (28) 

which was derived with the help of projection operator techniques from the same 
microscopic model using Born approximation in conjunction with a short memory 
approximation. As a main difference, in Agarwal's equation the K s term is absent. 
Here k is a phenomenological damping coefficient. 

From Eqs. (27) and (28) we obtain the following conditions for purity of states: 

K c a qq {t) + ^-* M (f) " \lc = (29) 



and, respectively, 



Mu coth(^)a qq (t) = ^. (30) 



The corresponding Schrodinger-type equations for a pure state have the Hamiltonian 

p 2 M 1 
H ' = 2M + + UJols ^ q2 + 2 lc(y<1P + ap ^ Q ~ ^(^P) ~ iK M - °"?(0) 2 

" Ss (p-a p (t))(q-a g (t)) (31) 



Muj 
and, respectively, 

H ' = 2M + + K(qP + apit)q ~ aq{t)p) ~ tKMuJ ° C ° th( 2^ )(g " aqit))2 - (32) 

4) All the above presented time-independent Liouville operators are not of Lind- 
blad form. In Ref. [17] it is shown that in the weak coupling limit, further coarse 
graining will result in a Lindblad operator. Indeed, for weak damping, the master 
equation (27) simplifies and takes on the following form, written in terms of usual 
creation and annihilation operators a), a: 



7. 

-i(u -r 

where 



p(t) = -t(uj + ^)[<fa,p(t)] + 7T ([atp(t), a ] + [a\p(t)a]) + 7i (M*),a t ] + [a, p(t)a%(33) 



^ = |MM|^)±1]. (34) 
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This equation, first derived by Weidlich and Haake [41] from a microscopic model for 
the damped motion of a single mode of the electromagnetic field in a cavity, is of 
Lindblad form and can be obtained formally as a particular case of the general master 
equation (68) for the damped harmonic oscillator (see next Section), if we take 

D pp = T ^Y 1 (li + 7t)> D m = ^J^ 7i + Dpq = °' X = ^ ll ~ 7T ^' ^ = °'( 35 ) 
From Eq. (33) we obtain the following condition for purity of states: 

2coth(^Kt a (t) = l (36) 
or, in terms of coordinate and momentum, 

(MwoMi) + ^ )coth( _^ H) , (37 ) 

For a pure state, the Schrodinger-type equation corresponding to Eq. (33) has the 
Hamiltonian 

H' = H + ih^(a a ,a - a J + 1) - ih^ C oth(^| ; )[(a t - a at )(a - a a ) + (38) 

with the notation H = Ti(ujq + ~f s /2)a)a. Taking into account the condition (37), we see 
that the mean values of the two Hamiltonians H and H' are equal: < H >=< H' > . 

In general, the dissipative systems cannot be described by pure states or by 
Schrodinger equations, because the environment produces transitions in any state 
basis. Nevertheless, we will show that this can happen in very limiting cases, cor- 
responding to certain special states. In order to find in the general Karrlein-Grabert 
model the states which remain pure during the evolution of the system, we consider the 
equations of motion for the second order moments of coordinate and momentum. To 
obtain these equations we first derive the evolution equation (17) with the coefficients 
(23) in coordinate representation: 

dp h 2 d 2 d 2 M lq 2 2 ih-f p d d 

+MD q (x -V)(fa + - \M 2 D v [x - yfp. (39) 

The first two terms on the right-hand side of this equation generate purely unitary 
evolution (with a renormalized potential). The third term is the dissipative term 
and the forth is the so-called "anomalous diffusion" term, which generates a second 
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derivative term in the phase space representation of the evolution equation, just like 
the ordinary diffusion term. The last term is the diffusive term, which is responsible 
for the process of decoherence. Since the considered dynamics is quadratic, we consider 
a density matrix solution of Eq. (39) of the form 

which is the general form of Gaussian density matrices (with zero expectation values of 
coordinate and momentum). By direct substitution of p into Eq. (39), we obtain the 
following system of equations satisfied by dispersions of coordinate and momentum: 



dcr qq (t) _ 2_ 
dt ~ .\/"< 



( 41 ) 



= -2 lp a pp (t) - 2M lq a pq (t) + 2M 2 D P , (42) 



da "" {i 1 -M lq a qq {t) + ^ pp (t) - lp a pq {t) + MD q . (43) 



dt " ww M 

Introducing the notation 



X(t) 



<rn&)/ m V% | ( 44 ) 

V <Tpq(t) 



and solving this system of equations with the method used in Refs. [6, 9], we obtain 
the solution: 

X(t) = T(X(0)-X(oo)) + X(oo), (45) 

where the matrix T is 

'fell fei2 fel 

T = | 621 b 22 fe 23 I , (46) 



e -7p< / °H °12 °13 
&31 ^32 ^33 



with time-dependent oscillating functions bij (i,j=l,2,3) given by (fl 2 = A^ q — 7^- 



- cos tit - 7 P — sin fit - 7,, (47) 
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bi2 = 7 9 (cosfit - 1), 



(48) 



bn = \/lq{lp cos m — fi sin fit — 7 P ), (49) 

621 =7?(cosfi*-l), (50) 

7 2 fi 
&22 = (y -7^)cosfit + 7 P — sinfit-7g, (51) 

^23 = \fTq{lp cos fit + fi sin fit - 7 p ) , (52) 

631 = -^(7? cos fit - I sin fit - (53) 

6 3 2 = -V^(y cos nt + § sin ^ - y )> ( 54 ) 

7 2 

& 33 = -27 g cosfit + -J. (55) 

The asymptotic values of variances and covariance have the following expressions: 

D p + lp D y M 2 D P 
<J qq {oo) = , a pp (oo) = , Opg(oo) = (56) 

7p7<? 7p 

a qq (0), (7pp(0), a pq (0). Introducing the expressions (23) for the coefficients D q and D p , 
we obtain the following equilibrium asymptotic values of the dispersions: 

<r qq (oo) =< q 2 >, o-pp(oo) =< p 2 >, a pq (oo) = 0. (57) 

If the asymptotic state is a pure state, then 

h 2 

a qq (oo)a pp (oo) =< q 2 >< p 2 >= — , (58) 

i. e. it is a squeezed state. Introducing the expressions of cr qq (t) and cr pq (t) given by 
(45) in the condition for purity of states (19), (23), we have shown, after a long, but 
straightforward calculation, that this condition is fulfilled, for any time t, only if the 
initial values of dispersions are equal to their asymptotic values: 

a qq (0) = a qq (oo), a pp (0) = a pp (oo), a pq (0) = a pq (oo). (59) 
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Then it follows from (45) that X(t) = X(oo), that is the dispersions have constant 
values in time, given by (57). Therefore, the state which preserves its purity in time is 
given by the density matrix (40), i. e. it is a squeezed state, with the well-determined 
constant variances a qq ,a pp (57). The fluctuation energy has also a constant value in 
time 

1 2 2 / x 

E= 2M <P >+V L< « >■ (60) 

At the same time, the total energy of the open system is given by the mean value of 
the Hamiltonian (16): 

< #o >= 2M app ® + ~2^ a ^ + + (61) 

and, since the expectation values of coordinate and momentum decay exponentially in 
time [19], the energy is dissipated and reaches the minimum value E. In the particular 
case of Agarwal model, the purity condition (30) shows that the variance of coordinate 
must also be constant in time: 

= H huM ' (62) 

2MM,coth(— ) 

Using this condition, we find from the equations of motion (41) - (43) written for the 
Agarwal model, when we have to take 

p = ~W ^2k^T^ q ' lp = K ' lq = ^ ' 



that the dispersions have to satisfy the following equalities: 

M*) = ^^coth(^), a pp (t) = M 2 u 2 a qq (t), a pq (t) = 0. (64) 



The relations (62), (64) are compatible only if coth.(huo/2kBT) = 1, that is only when 
the temperature of the thermal bath is T = 0. Then finally we get 

h hMu 
^ = 2AW <7 «' = -2 _ ' <7 " = (65) 

and, therefore, in the particular case of Agarwal model, the usual coherent state is 
the only state which remains pure for all times, if the temperature is T = 0. In this 
case the fluctuation energy of the harmonic oscillator keeps all the time its minimum 
value E min = Kuoq/2. The relation (61) shows that in this case the total energy is 
also dissipated and reaches E min . The same results can be obtained for the model of 
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Weidlich and Haake, described by the evolution equation (33). Indeed, this model 
is a particular case (cf. Eqs. (35)) of the Lindblad model considered in the next 
Section and from the purity condition (37) it follows again that the coherent state is 
the only state which preserves its purity during the evolution in time of the system, for 
a zero temperature of the thermal bath. The importance of the states which preserve 
their purity in time will become evident in Sec. 5, in the context of discussing the 
decoherence phenomenon. 



4 Purity of states in the Lindblad model 

We now consider the condition for purity of states in the Lindblad model for the 
damped harmonic oscillator, based on quantum dynamical semigroups [2, 3, 7, 8]. 
The most general Markovian evolution equation preserving the positivity, hermiticity 
and trace of p can be written as: 



dp(t) 



dt ~[H, P (t)} + ^([V jP (t),V}} + [% P {t)V}\). (66) 

j 

Here H is the Hamiltonian operator of the system and Vj, Vj are operators on the 
Hilbert space Ti. of the Hamiltonian which model the interaction with the environment. 
In the case of an exactly solvable model for the damped harmonic oscillator we take the 
two possible operators V\ and Vi linear in p and q [6, 8, 9] and the harmonic oscillator 
Hamiltonian H is chosen of the general form 

H = Ho + \ (qp #o = + ^fq\ (67) 

With these choices the master equation (66) takes the following form [6, 9]: 

dp i * i % 

-j t = -^[H ,p\ - —(A + p)[q,pp + pp] + ^(A-^)[p,pg + gp] 

-fm m - m + fFa?, m + % %p\\\ (68) 

The quantum diffusion coefficients D pp , D qq , D pq and the dissipation constant A satisfy 
the following fundamental constraints [6, 9]: D pp > 0,D qq > and 

D PP D qq - Dl q > ^f. (69) 

The relation (69) is a necessary condition that the generalized uncertainty inequality 
(2) is fulfilled. 
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By using the complete positivity property it was shown in [9] that the relation 
Tr(p(t) £ VjVj) = ]T Tr(p(t)Vj)Tr(p(t)Vj) (70) 

3 3 

represents the necessary and sufficient condition for p(t) to be a pure state for all 
times t > 0. This equality is a generalization of the pure state condition [42-44] to 
all Markovian master equations (66). If p 2 {t) = p(t) for all t > 0, then there exists 
a wave function ip G H which satisfies a nonlinear Schrodinger equation with the 
non-Hermitian Hamiltonian 

H f =H+ij2< m, v}m >vj- l -< m,n v^m >~y. v}% ( 71 ) 

3 3 3 

For the damped harmonic oscillator the pure state condition (70) takes the form [9] 

h 2 \ 

D. pp cr qq {t) + D qq a P p(t) - 2D pq a pq {t) = (72) 

and the Hamiltonian (71) becomes 

H' = H + \(a p (t)q - a q (t)p) - ^[D pp (q - a q (t)) 2 + D qq {p - a p {t)) 2 

-D pq ((p - a p {t)){q - a q {t)) + (q - <j q (t))(p - a p (t))) - —]. (73) 

It is interesting to remark that the mean value of this Hamiltonian in the state p(t) 
is equal to the mean value of the Hamiltonian H. From a physical point of view this 
result is quite natural, since the average value of the new Hamiltonian H' describing 
the open system must give the energy of the open system. 

In order to find the Gaussian states which remain pure during the evolution of the 
system for all times t, we consider the pure state condition (72) and the generalized 
uncertainty relation for pure states: 

h 2 

<rpp(t)a qq (t) - a 2 pq (t) = -. (74) 
By eliminating a pp between the equalities (72) and (74), like in [45], we obtain: 

/ / 1\ DpqVpqjt) + jfc 2 \ 2 DppDgg ~ ^pq w /.x J^XDpq ^ 2 

K a qq\ l ) n ) "I - n2 l{°pqy l ) n n _ r)2 > 

J-Spp J^pp ^ppJ^qq J-Spq 

1 n 2 D PP D m - D lg- l k2X2 D n l_o m\ 

+ A (n n _ D 2\2 UppL) qq \-\). [(b) 
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Since the diffusion and dissipation coefficients satisfy the inequality (69), we obtain 
from Eq. (75) the following relations which have to be fulfilled at any moment of time: 

DppD qq - — ^-p-, (76) 



h 2 \ 

Dp P o qq {t) - D pq a pq (t) — = 0, (77) 



a pq (t)(D pp D qq - D 2 pq ) - ^D pq = 0. (78) 

From relations (74) and (76) - (78) it follows that the pure states remain pure for all 
times only if the variances have the form: 

M*) = ^ = *«(*) = (79) 

i. e. they do not depend on time. If these relations are fulfilled, then the equalities 
(72), (74) and (76) are equivalent. Using the asymptotic values of variances for an 
underdamped oscillator (given by Eqs. (3.53) in [9]) and the relations (79), we obtain 
the following expressions of the diffusion coefficients which assure that the initial pure 
states remain pure for any t (f2 2 = uo 2 — /i 2 ): 

D « = ^ D » = ^T> D ™ = -^T- (80) 

Formulas (80) are generalized Einstein relations and represent typical examples of 
quantum fluctuation-dissipation relations, connecting the diffusion with both Planck's 
constant and damping constant [4, 46]. With (80), the variances (79) become 

h Timu 2 U\x 

Then the corresponding state described by a Gaussian Wigner function is a pure 
quantum state, namely a correlated coherent state [18] (squeezed coherent state) with 
the correlation coefficient (4) r = —/i/cu. Given a qq , a pp and <r pq , there exists one 
and only one such a state minimizing the uncertainty a (2) [47]. A particular case of 
Lindblad model (corresponding to A = /i and D pq = 0) was considered by Halliwell and 
Zoupas by using the quantum state diffusion method [22] . We have considered general 
coefficients A and /i and in this respect our expressions for the diffusion coefficients 
and variances generalize also the ones obtained by Dekker and Valsakumar [45] and 
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Dodonov and Man'ko [48], who used models where X — /i was chosen. If /i = 0, we 
get D pq = from (80). This case, which was considered in [10], where we obtained 
a density operator describing a pure state for any t, is also a particular case of our 
present results. For /i = 0, the expressions (81) become 

^ ft muJ n /qo\ 

a qq = °vv = > a pg = °» ( 82 ) 

which are the variances of the ground state of the harmonic oscillator and the corre- 
lation coefficient is r = 0, corresponding to the usual coherent state. 
The fluctuation energy of the open harmonic oscillator is 

E{t) = ^ PP (t) + \mu 2 (J qq (t) + fJUTpgit). (83) 

If the state remains pure in time, then the variances are given by (79) and the fluctu- 
ation energy is also constant in time: 

E = \(^ D pp + \rruv 2 D qq + fjiDpq). (84) 

Minimizing this expression with the condition (76), we obtain just the diffusion coef- 
ficients (80) and E min = hil/2. Therefore, the conservation of purity of state implies 
that the fluctuation energy of the system has all the time the minimum possible value 
E m i n . The total energy of the open system is given by the mean value of Hamiltonian 
(67): 

1 ~o fnu 2 ~ u 

< H >= < p z > H < q z > +- < qp + pq > 

2m 2 2 

= 2^°PP(*) + 2 mUj2(J 1l( t ) + ^pqif) + + ~^ mW ' 2(T W t ) + VVp(t)(J q (t) (85) 

and, since the expectation values of coordinate and momentum decay exponentially in 
time [6, 9], the energy is dissipated and reaches its minimum value E min . 

If the asymptotic state is a Gibbs state [6, 9], then the condition (76) on the 
diffusion coefficients is satisfied only if /i — and the temperature of the thermal 
bath is T = 0. Like in the Agarwal and Weidlich-Haake models, discussed in the 
previous Section, in this limiting case the influence on the oscillator is minimal and 
Emin = hw/2, which is the oscillator ground state energy, the correlation coefficient (4) 
vanishes and therefore the correlated coherent state (squeezed coherent state) becomes 
the usual coherent (ground) state. 

The Lindblad equation with the diffusion coefficients (80) can be used only in 
the underdamped case, when u > ji. Indeed, for the coefficients (80) the fundamental 
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constraint (69) implies that m 2 {u 2 — ii 2 )D 2 qq > h A 2 /4, which is satisfied only if u > /i. 
It can be shown [48] that there exist diffusion coefficients which satisfy the condition 
(76) and make sense for ui < //, but in this overdamped case we have always a > h 2 /4 
and the state of the oscillator cannot be pure for any diffusion coefficients. 

If we choose the coefficients of the form (80), then the equation for the density 
operator can be represented in the form (66) with only one operator V, which up to a 
phase factor can be written in the form: 



i/ = ] jw; q [{ Y~ iDpq) « + iDqqp)l = 2RA - (86) 

The correlated coherent states (6) with nonvanishing momentum average can also 
be written in the form: 

= ( 2^^ eXp[_ 4^ (1 " f apq) ^ X ~~ ^ + Ji apX] (87) 
and the most general form of Gaussian density matrices compatible with the general- 
ized uncertainty relation (2) is the following: 

,., .1,1 r 1. , X ~\~ 1J / \ \ 9 

Zncr qq la qq I 

2 

+ ^(^ " ° q ){x -y)- ^(a pp yf + ^(x - y)\. (88) 

These matrices correspond to the correlated coherent states (87) if a qq , a pp and a pq in 
(88) satisfy the equality (2), in particular if the variances are taken of the form (81). 
Consider now the harmonic oscillator initially in a correlated coherent state of the 
form (87), with the corresponding Wigner function (7). For an environment described 
by the diffusion coefficients (80), the Wigner function at time t is given by 

W(p,q,t) = ^- 
7m 

2 

x exp{-^ [cr pp (g - (T q (t)) 2 + cr qq (p - a p (t)) 2 - 2a pq (q - cr q (t))(p - a p (t))]}, (89) 

with the constant variances (81). The correlated coherent state (squeezed coherent 
state) remains a correlated coherent state with variances constant in time and with 
<J q (t) and a pit) giving the average time-dependent location of the system along its 
trajectory in phase space. In the long-time limit <r q (t) = 0, <r p (t) = and then we have 

< x\p(oc)\y >= (^)^xp{-g[fi(x 2 + y 2 ) + i^x 2 - y 2 )]}. (90) 
The corresponding Wigner function has the form 

HWP, ,) = -exp[--(- + -.W + m)\- (91) 
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5 Entropy and decoherence 



Besides the von Neumann entropy S (13), (14), there is another quantity which can 
measure the degree of mixing or purity of quantum states. It is the linear entropy Si 
defined as 

^ = Tr(p-p 2 ) = l-Trp 2 . (92) 

For pure states Si — and for a statistical mixture Si > 0. As it is well-known, the 
increasing of the linear entropy Si (as well as of von Neumann entropy S) due to the 
interaction with the environment is associated with the decoherence phenomenon (loss 
of quantum coherence), given by the diffusion process [23, 24]. Dissipation increases 
the entropy and the pure states are converted into mixed states. The rate of entropy 
production is given by 

Si(t) = -2Tr(p'p) = -2Tr(p£(p)), (93) 

where C is the evolution operator. According to Zurek's theory [23, 24], the maximally 
predictive states are the pure states which minimize the entropy production in time. 
These states remain least affected by the openness of the system and form a "preferred 
set of states" in the Hilbert space of the system, known as the "pointer basis". Deco- 
herence is the mechanism which selects these preferred states - the most stable ones 
under the evolution in the presence of the environment. 

For the models of the damped harmonic oscillator considered in this paper, we 
can obtain the expressions for the rate of entropy production given by Eq. (93). For 
Gaussian states the linear entropy (92) becomes 

S^) = l-I v=\y/d (94) 
and then the time derivative of the linear entropy is given by 

1 dv _ h da qq {t) da pp (t) 9 ^WX> mi (Qt * 

Si{t) - ^ - — [—^-*„{t) + —^-<r„{t) - 2—a„{t)]. (95) 

From the system of equations (41) - (43) for the Karrlein-Grabert model we obtain 

Si{t) = ^=[M 2 D p (t)a qq (t) - MD q it)a pq it) - lp {t)a}. (96) 

Suppose at the initial moment of time t — the state is pure. When the conditions 
(19), (20) for purity of states are fulfilled for all t, the expression of the rate of linear 
entropy becomes 

Siit) = ^[M 2 D p {t)a qq {t) - MD q {t)a pq {t) - ^ 7p (i)] = (97) 
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and then the entropy production has its minimum value Si — 0. For the thermal initial 
condition with the coefficients (18), the rate of entropy production is given by 

^ l{t) = 2^ [lp{t) <p2> aqq{t) ~ {Mlq{t) <q2> ~^r )apq{t) ~ 7p(t)a] ' (98) 

for strictly Ohmic damping it is 

^ = 2^ [l <p2> aqq{t) ~ {Muj2 ° <q2> ~^r )apq{t) ~ la] (99) 

and for Drude damping the rate of entropy production is also given by an expression 
like (98), where now 7 P = 2a and ^ q = a 2 +r] 2 . When the condition for purity is fulfilled 
for any t, the values of the rate of linear entropy given by (98), (99) become also 0. 
According to the results of Sec. 3, if the condition for purity of states is fulfilled for 
any t in the Karrlein-Grabert model, then the Gaussian state will be a pure squeezed 
state, with constant in time variances. At the same time the rate of linear entropy 
production vanishes and, therefore, according to the Zurek's theory of decoherence, 
the most stable states are the pure squeezed states, with constant variances. The same 
conclusion is valid for the weak damping model, given by the master equation (27), 
for which the rate of entropy production has the expression 

s<« = 5S [ ^ M ' ) + ^ M(t) ~ 7 ^ (100) 

while for the Agarwal model given by the master equation (28) we obtain 

Analogously, for Eq. (33) of Weidlich and Haake, the rate of entropy production is 
given by 

AW ^ [( M.„M f ) + ^)coth(^)-f] (102) 

and, according to the results of Sec. 3, for Agarwal and Weidlich-Haake models, the 
usual coherent states are the most stable ones under evolution in the presence of the 
environment. Using Eq. (93) for the Lindblad equation (68), we obtain the following 
rate of entropy production: 

S t (t) = ^[D pp Tr(p 2 q 2 - pqpq) 

h 2 \ 

+D qq Tr(p 2 p 2 - pppp) - D pq Ti(p 2 (qp + pq) - 2pqpp) - — Tr(p 2 )] (103) 
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[D pp a qq (t) + D qq a pp (t) — 2D pq a pq (t) — 2\a\. 



(104) 



If the initial state is pure, then according to the complete positivity property of the 
Lindblad model we have 



which means that the linear entropy can only increase, so that the initial pure state 
becomes mixed. When the state remains pure, Eq. (104) becomes, cf. Eq. (72) : 



and, therefore, the entropy production will be Si = 0. Since the only initial states 
which remain pure for any t are the correlated coherent states, we can state that in the 
Lindblad theory these states are the maximally predictive states. The present results, 
obtained in the framework of Karrlein-Grabert and Lindblad models, generalize the 
previous results which assert that for many models of quantum Brownian motion in 
the high temperature limit the usual coherent states correspond to minimal entropy 
production and, therefore, they are the maximally predictive states. As we have seen, 
such coherent states can be obtained in the Lindblad model as a particular case of 
the correlated coherent states by taking /i = 0, so that the correlation coefficient (4) 
r = 0. Namely, Paz, Habib and Zurek [23, 24] considered the harmonic oscillator 
undergoing quantum Brownian motion in the Caldeira-Leggett model and concluded 
that the minimizing states which are the initial states generating the least amount 
of von Neumann or linear entropy and, therefore, the most predictable or stable ones 
under evolution in the presence of an environment, are the ordinary coherent states. 
Using an information-theoretic measure of uncertainty for quantum systems, Anderson 
and Halliwell showed in [25] that the minimizing states are certain general Gaussian 
states. Anastopoulos and Halliwell [26] offered an alternative characterization of these 
states by noting that they minimize the generalized uncertainty relation. According 
to this assertion, we can say that in the Lindblad model the correlated coherent states 
are the most stable ones which minimize the generalized uncertainty relation (2). 
Our result confirms that one of [26], where the model for the open quantum system 
consists of a particle moving in a harmonic oscillator potential and linearly coupled 
to an environment consisting of a bath of harmonic oscillators in a thermal state. We 
remind that the Caldeira-Leggett model considered in [23, 24] violates the positivity 



5,(0) = ^[D pp a qq (0) + D qq a pp (0) - 2D pq a pq (0) - — ] > 0, 



(105) 



Si(t) = j^[D pp a qq (t) + D qq a pp (t) - 2D pq a pq (t) — ] = 



(106) 
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of the density operator at short time scales [49, 50], whereas in the Lindblad model 
the property of positivity is always fulfilled. 

The rate of predictability loss, measured by the rate of linear entropy increase, is 
also calculated in the framework of Lindblad theory for the damped harmonic oscillator 
by Paraoanu and Scutaru [27], who have shown that, in general, the pure or mixed 
state which produces the minimum rate of increase in the area occupied by the system 
in the phase space is a quasi-free state which has the same symmetry as that induced 
by the diffusion coefficients. For isotropic phase space diffusion, coherent states (or 
mixture of coherent states) are selected as the most stable ones. In order to generalize 
the results of Zurek and collaborators, the entropy production was also considered by 
Gallis [28] within the Lindblad theory of open quantum systems, treating environment 
effects perturbatively. Gallis considered the particular case with D pq = and found 
out that the squeezed states emerge as the most stable states for intermediate times 
compared to the dynamical time scales. The amount of squeezing decreases with time, 
so that the coherent states are most stable for large time scales. For D pq ^ our 
results generalize the result of Gallis and establish that the correlated coherent states 
are the most stable ones under the evolution in the presence of the environment. 

6 Summary 

In the present paper we have first considered the generalized quantum master equa- 
tions derived by Karrlein and Grabert [17] for the microscopic model of a harmonic 
oscillator coupled to a harmonic bath. We have obtained the conditions for purity of 
states for different initial conditions and different types of damping, including strictly 
Ohmic, Drude and weak coupling cases, Agarwal and Weidlich-Haake models. We have 
shown that the states which remain pure all the time are the pure squeezed states with 
well-determined constant in time variances. For pure states, we have also derived the 
corresponding generalized Schrodinger-type nonlinear equations. Then we have stud- 
ied the one-dimensional harmonic oscillator with dissipation within the framework of 
Lindblad theory and have shown that the only states which stay pure during the evo- 
lution of the system are the correlated coherent states, under the condition of a special 
choice of the environment coefficients, so that the variances and covariance are con- 
stant in time. We have also obtained the expressions for the rate of entropy production 
in the considered models and have shown that the states which preserve their purity 
in time are also the states which minimize the entropy production and, therefore, they 
are connected with the decoherence phenomenon. According to the Zurek's theory 
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of decoherence, in Karrlein-Grabert and Lindblad models, as well as in the consid- 
ered particular models, these states are the most stable ones under the evolution of 
the system in the presence of the environment. In a next work in the framework of 
these theories we plan to discuss in more details the connection between uncertainty, 
decoherence and correlations of open quantum systems with their environment. 
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